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THE  USE  OF  VECTOR  FIELDS  TO  MODEL  THE  PHYSICAL  BLOCKAGE 
FROM  POWER  SUPPLY,  CABLE,  AND  TRANSFORMER  SOURCES 


INTRODUCTION 

This  effort  was  designed  to  improve  the  electromagnetic  (EM)  analysis  techniques 
that  can  be  used  to  supplement  the  Intelligent  EMC  Analysis  and  Design  System 
(lEMCADS)l  developed  at  the  Electromagnetic  Compatibility  (EMC)  Branch  of  the 
Naval  Undersea  Warfare  Center  Detachment  in  New  London,  Connecticut.  The  primary 
objective  of  this  work  was  to  model  the  physical  blockage  of  magnetic  fields  produced  by 
the  lEMCADS  sources. 

The  lEMCADS  uses  such  sources  as  transformers,  cables,  and  power  supplies, 
which  it  represents  in  a  simplistic  and  scalar  form  to  yield  worst-case  magnetic  field 
conditions.  These  sources  do  not  give  rise  to  uniform  fields  or  plane  waves,  but  to  fields 
that  vary  as  inverse  powers  of  the  radial  distance. 

It  was  decided  that  the  best  way  to  obtain  an  analytical  solution  to  this  problem 
was  to  expand  the  sources  that  exist  in  lEMCADS  in  Legendre  polynomials  and  use  a 
hollow  cylindrical  physical  blockage  as  a  representation  of  cabinets.  This  technique 
enables  the  fields  in  the  various  regions  around  and  within  the  physical  blockage  to  be 
modeled  as  expansions  in  Legendre  polynomials. 2-6  These  expansions  make  the  source 
modeling  compatible  with  the  physical  blockage  modeling,  and  thus  boundary  conditions 
can  be  applied  for  the  B  and  H  fields*  at  the  two  interfaces  of  the  cylindrical  shell.  The 
application  of  the  boundary  conditions  generates  a  system  of  10  equations  and  10 
uiiknowns  for  the  coefficients  of  expansion  for  each  source  type.  The  solution  of  the  10 
equations  and  10  unknowns  yields  the  coefficients  of  expansion  of  the  Legendre 
polynomials.  The  final  Legendre  polynomial  expansions  produce  the  magnetic  scalar 
potentials  for  the  different  sources.  The  gradients  of  the  magnetic  scalar  potentials  are 
then  taken  to  yield  the  vector  B  and  H  fields  in  the  regions  of  the  source,  the  regions 
outside  and  within  the  shell  of  the  physical  blockage,  and  the  hollow  region  of  the 
physical  blockage. 


*  According  to  Feynman's  definition,  B  represents  the  magnetic  field  vector  and  H  represents  the 
magnetizing  field. 
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Point  of  Observation 
(Measurement) 


THE  MATHEMATICAL  MODELS 


The  lEMCADS  power  supply  field  varies  as^^_p/|  and  the  lEMCADS 
transformer  varies  as  _  -,|3 .  From  figure  1,  r'  is  the  distance  from  the  source  to  the 
center  of  the  cylindrical  shell  and  f  is  the  distance  from  the  cylindrical  axis  to  the  point 
of  observation.  The  coordinates  of  the  point  of  observation  (experimental  point)  are 
always  rand  6. 


GENERAL  METHOD  FOR  APPLICATION  OF  THE 
BOUNDARY  CONDITIONS 

The  coefficients  of  the  Legendre  polynomial  expansions^  are  solved  by 
application  of  the  boundary  conditions  at  the  two  interfaces:  r  =  a  and  r  =  Zj.  These 
boundary  conditions  are  the  continuity  of  the  tangential  components  of  the  H  field  and 
the  normal  component  of  the  B  field.  To  apply  these  boundary  conditions  requires  that 
the  coefficients  of  ^cosine^^  j^sine°^j,  where  n  =  0,1,2  and  m  =  0,1,2,  be  set  equal  to  zero. 

The  application  of  the  boundary  conditions  involves  two  cases:  r>f'  and 
f  <  f '.  The  coefficients  of  expansion  of  the  Legendre  polynomials  are  different  for  these 
two  cases  since  the  sources  are  expanded  as  a  binomial  expansion  in  powers  of  for 
the  case  r>r'and  in  powers  of  for  the  case  r<r' .  (Physically,  the  case  r>r' 
means  that  the  source  is  closer  to  the  physical  blockage  than  the  point  of  observation  and 
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the  case  r'>r  means  that  the  point  of  observation  is  closer  to  the  physical  blockage  than 
the  source.)  The  binomial  expansions  of  the  sources  are  necessary  to  obtain  coefficients 
of  powers  of  sines  and  cosines.  These  expansions  can  then  be  combined  with  the 
Legendre  polynomial  expansions  to  generate  the  10  equations  and  10  unknowns.  The  10 
equations  yield  unique  solutions  for  the  10  unknowns,  which  are  then  substituted  back 
into  equations  (1),  (2),  and  (3)  (see  below).  The  gradient  operations  were  carried  out  in 
these  three  equations  to  yield  the  general  solutions. 

THE  MATHEMATICAL  FORM  FOR  THE  }\r-r\  SOURCE 

The  H  fields  are  given  in  three  regions:  r>b,  a<r<b,  and  r<a.  For  r  >  b,  the 
H  field  is  given  by 


g  (r-r) 


r-r  r-r 


^=0 


(1) 


For  a<r<b. 


Pf  (cos(0))J 


(2) 


and  for  r<a, 

H  =  s/p,  {com)-  (3) 

^=0 


THE^  _  SOURCE  IN  COMPONENT  FORM  AND  ITS  BINOMIAL 
EXPANSION 

The  B  field  for  the  can  be  expressed  as 
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where  the  second  factor  is  the  unit  vector  extending  from  the  source  to  the  observation 
point.  Taking  the  dot  product  of  B  with  the  unit  vectors  in  the  radial  and  angular 
directions  of  the  physical  blockage  cylindrical  coordinate  system  yields 


B  =  — ^•-=  ^-~^{r-r'cos6) 

r  s-  s'\  s-  S-'  =•  2  V  / 


\r  -r  \  \r  -  r  \  r  \r  -r 


and 


B  ^  rs  r'smd 

r -r  I  r-r 


Expanding  _,|2  in  the  binomial  expansion  yields 


-  -/2  ~ 
r-r  r 


^  r'^  2r' 

1 - ^  + - COS0 

r  r 


for  f>f' 


and 


r  2r 
\  ^  ^  J 


for  r<r' 


(4) 


(5) 


(6) 


(7) 


THE  MATHEMATICAL  SOLUTION  OF  THE  B  AND  H  FIELDS  FOR  THE 
SOURCE 

We  substitute  equations  (6)  and  (7)  into  equations  (4)  and  (5)  and  the  result  into 
the  source  terms  of  equations  (1),  (2),  and  (3).  By  applying  the  boundary  conditions  for 
the  B  and  H  fields  to  the  resulting  equations  and  setting  the  coefficients  of  the  powers  of 
the  sines  and  cosines  equal  to  zero,  we  generate  the  10  equations  and  10  unknowns  and 
solve  them  using  the  Mathematica  7  software,  as  shown  in  the  appendix  of  this  report. 

The  boundary  aXr  =  b  has  two  cases:  f<f'  and  f  >  f ' .  For  the  continuity  of 
B^,  the  coefficients  of  {cos  6)^  are 
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(8) 


A. 


r  u^\ 
b  — 


.'2 

V  '  J 


1  3  1 


1 


A,, 


=  +  for  r<f' 


and 


/2\ 


4U-- 

Z) 


13  1  1  o  o  I.  3  7,  , 

for  r>r'; 


(9) 


the  coefficients  of  (cos0)^  are 


r' 


-r 
V  '  J 


=  _2ai^-2^7i^  for  r  <r 


(10) 


and 


r^>^  ^ 

T  +  r' 

y 


=  -2«iiJ-2/^ri7j  for  r>f'; 

D  0 


(11) 


and  the  coefficients  of  (cos  6)  are 


=  for  r<r’  (12) 

and 

-Ir'^B  9  1  9  r 

— 7r^  =  T«2— +  3/i^2^--^^  for  r>f'.  (13) 

Dir  Zb 

For  the  continuity  of  Hq  ,  there  are  no  terais  in  the  coefficient  for  (cos  0)®.  For  the 
coefficient  of  (sin)^. 


f•  =  -A-f  for  r-<f' 


(14) 


and 


A 

b'^ 


^  =  -A-7i-  for  F>r'. 


(15) 
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For  the  coefficient  of  (sin  B  cos  0)^, 


^^0^  ,  3ot2  u2  1  ^7; 


,2  =  for  f<f' 


and 


IBy^  Sa.  0/3  1.2  37,  . 

_^  +  _^  =  3/32h"+-^  for  r. 


For  the  continuity  of  at  r  =  a,  the  coefficient  of  (cos0)°  is 


U7o  o  37,  d.a 
2  r-^2  2  2 


the  coefficient  of  (cos  6)^  is 

and  the  coefficient  of  (cos  0)^  is 

|if3^2^-|4]  =  |4a  • 

V  2  a  J  2 


For  the  continuity  of  H  at  r  =  a,  the  coefficient  of  (sin)i  is 


Aa+7i  — =  V 


a 


and  the  coefficient  of  (sin  B  cos  0  )Hs 


^2^  "^72  3  ^2^  • 


a 


The  radii  a  and  h  are  fixed  known  numbers  and  the 
ocq,  cci,a2, Pi, P2,  To >71^  72’^!’  ^2’  tho  following  solutions: 


(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

unknowns  are 
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+18/iflV  +12;i  Vz?^  +4/i6^  +  6/iVj 

HV  V 


6X2X4^ 


-6b'^X2X^  -  SiJX^Xsr'  -  (t^^X^X^r'  -  ^bX^X.^/^  +  eb'^X^Xf^  r''^ ) 


6X2^4'- 


A 


«!  = 


_B„{-3b^  +  r'^) 


2r 


/3  » 


a,  = 


24^ 


^  3X2  r'2  ’ 


_  _B^{l-3pi  +  2\i^)a^b^X, 
|iX2r'^ 


S(l-^)^ 
2X,r'^  ’ 


72  = 


2B^{-l  +  2ll)a^b^X^ 

3X,r'^ 


^  2X2  r'^  ’  ^  3X2  r'2  ’ 


which  were  solved  numerically  using  the  Mathematica  software  (see  the  appendix). 


Substituting  the  above  Legendre  polynomial  coefficients  back  into  equations  (1), 
(2),  and  (3)  yields  for  r  >  h 


H  =  f\ 


^ - (r  -  r' cos  0)- -^0080+ -^+^^(3008^  0-l] 

.r  +r'  —  27t'  r  r  2r  ' 


+  0 


BrsAsind  Oi  .  ^  a,  ^ 

^  2 - 1-  sm  0  +  — X  (6  cos  0  sin  0) 


r^  +  r'^-2rr'  r 


2r 


for  b>r> a 


1 


H  =  r 


7  -  A  ® ^ 7  2^^r](3cos^  e  - 1) 


+  e\ 


(  y  ^ 

V  ^  J 


67 

sin  0  -  +  — ^  sin  0  COS  0 

z  r 


and  for  r  <  Z> 


1 

'2S  3S  ,  ^  ,  S' 

■5,  35, 

II 

7rcos6>  4^U3cos  e  1)+^ 

+  e 

-4  sin  0  +  — ^  sin  0  cos  6 
r  r 

THE  MATHEMATICAL  FORM  FOR  THE  ^  -  7'  SOURCE 

The  H  fields  are  given  in  three  regions:  r>h,  a<r<b,  and  r<a.  For  r>b,  the 
H  field  is  given  by 


H  = 


^  (r-r') 

\r-r'f  \f-r\ 


2 

vs- 


(X 


(cos(0)) ; 


for  a<r<b, 

S=-vxfj3/'+i^  Me))' 

i=oV  ^  J 

and  for  r<  a, 

h=-vY,5/p^  (cos(e)). 

1=0 


(23) 


(24) 


(25) 


Expanding  equations  (23),  (24),  and  (25)  yields  for  r  >  h 


H  =  f 


°  (r-r'cosO) — ^cos0  +  — ^  +  -^— ^(3cos^  0  — l) 


^^2  —  Orr']  ^ 


-2rr') 


+  d 


Br' sinO  a,  cc 

— - 2 — ^sin0  +  — 7(6cos0sin0) 

.-2  ,  ^/2  ^  2r 


{r^+r'^-2rr')  '' 
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for  b>r> a. 


H  =  r 


7  0  +  ^  +  (r2  7  +  2/?/ j(3cos^  d-l) 


+  e 

and  for  r  <b, 

\ 

H  =  f\ 


^  y  ^  67, 


P  +-^  \sm6-p  r  +  —^sin0cos0 


\  '  J 


^^cos0--^^(3cos^  0-1)+-^ 
r  Ar  r  . 


-I  5a  35o 

+  01  -jsin0  +  — ^sin0cos0 
r  r 


THE  /^r-rf  SOURCE  IN  COMPONENT  FORM  AND  ITS  BINOMIAL 
EXPANSION 


The  B  field  for  the 


can  be  expressed  as 


where  the  second  factor  is  the  unit  vector  extending  from  the  source  to  the  observation 
point.  Taking  the  dot  product  of  B  with  the  unit  vectors  in  the  radial  and  angular 
directions  of  the  physical  blockage  cylindrical  coordinate  system  yields 


-r'cos  0) 


and 


5  = 


r  -  r 


(26) 


(27) 


Expanding  -/|4  ^  binomial  expansion  yields 
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1  f  fr'^  2r'  1 

- for  r>r 

^\r-r'\  /  r 


X-  -.14^-J  1-2 -72-^cos0  for  r<r' 
\r-r\  /  \  [r'  r  ) 


THE  MATHEMATICAL  SOLUTION  OF  THE  B  AND  H  FIELDS  FOR  THE 
Xr-r\^  SOURCE 


We  substitute  equations  (28)  and  (29)  into  equations  (26)  and  (27)  and  the  result 
into  the  source  terms  of  equations  (23),  (24),  and  (25).  By  applying  the  boundary 
conditions  for  the  B  and  H  fields  to  the  resulting  equations  and  setting  the  coefficients  of 
the  powers  of  the  sines  and  cosines  equal  to  zero,  we  generate  the  10  equations  and  10 
unknowns,  which  are  solved  to  obtain  the  coefficients  of  expansion.  These  coefficients 
are  then  substituted  back  into  equations  (1),  (2),  and  (3)  to  yield  a  numerical  solution. 

The  boundary  air  =  b  has  two  cases:  f<r'  and  f>f'.  For  the  continuity  of  , 
the  coefficients  of  (cos  6)^  are 


■^0  1  _.r2_ 


2b^  ]  a.  Sf  a. 


+  for  r<f'  (30) 


bJ  2r'^]  ccq 

b^[  b^  j  b^  2W 


}\  =  -^yo\i2\-l^p2b  +  j[^\  for  r>r';  (31) 


the  coefficients  of  (cos0)^  are 


Bj4b^  ^  2b^]  2a,  „  ( \  r 
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E  (  1  2a  f  2\ 


and  the  coefficients  of  (cos  $)^  are 


=  for  Fsr 


■4(4'-")-^-  =  3ff2ft6-^l  for  F>r 


For  the  continuity  of  at  r  =  b,  there  are  no  terms  in  the  coefficient  for  (cos  Of.  For 
the  coefficients  of  (sin)i. 


r  \  • 


A^+7i^J  for  f<r 


rM 


=  -  A^  +  7iTT  for  r^r' 


For  the  coefficients  of  (sin  B  cos  0)^ 

=  for  r.r 


=  torr-.r. 


For  the  continuity  of  at  r  =  a,  the  coefficient  of  (cos  B'f  is 


2^2^  -  72  4-1 = 0 ; 


of-'  2 


the  coefficient  of  (cos  B^  is 


11 


and  the  coefficient  of  (cos  Of  is 


I- 1  =  35,0, 
2  V 


For  the  continuity  of  at  r  =  a,  the  coefficient  of  (sin)^  is 
-(j5,a  +  r,j^^  =  -5,a 
and  the  coefficient  of  (sin  9  cos  0)Hs 


(42) 


(43) 


(44) 


THE  MATHEMATICAL  SOLUTION  OF  THE  B  AND  H  FIELDS  FOR  THE 

y^r-r\  cable  source 


The  cable  source  is  perpendicular  to  the  radial  field  source  and  is  of  the  same 
magnitude.  The  cable  source  field  in  the  presence  of  the  physical  blockage  was  solved 
with  the  binomial  and  Legendre  polynomial  expansions  in  the  same  way  that  the 
solutions  for  the  previous  two  sources  presented  in  this  report  were  obtained.  The 
solution  to  the  fields  at  0  =  0  yields  disturbance  fields,  which  are  functions  of  angular 
derivatives  of  the  Legendre  polynomials.  The  angular  derivatives  of  the  Legendre 
polynomials  yield  the  angular  component  of  the  disturbance  field,  and  each  term  has  a 
sin  6  factor.  Setting  6=0  yields  zero  for  the  disturbance  field,  which  the  experiments 
conducted  at  0  =  0  proved  to  be  true,  as  shown  in  figure  2. 


The  source  field  from  the  cable  is  given  by 


Bw=(  %-r’\ 


where  8^=  —  (I  is  the  current  in  the  wire)  and  6^  is  the  unit  vector  in  the  angular 


direction. 


Ik 


It  can  easily  be  seen  that  the  radial  field  is  perpendicular  to  the  angular  field  so 

that 

Rewriting  B^  in  the  coordinate  system  of  the  physical  blockage  yields 


K  =  ’ 

and  rewriting  B  of  the  radial  field  in  the  coordinate  system  of  the  physical  blockage 
yields 

B  =  B,f  +  Bse  . 

Taking  the  dot  product  of  the  two  previous  equations  results  in 

Since  the  sources  are  orthogonal  fields,  the  fields  around  the  physical  blockage 
are  orthogonal  as  well.  (This  property  is  also  characteristic  of  potential  flow  where 
streamlines  (velocity  fields)  and  equipotential  surfaces  are  orthogonal  and  the  Cauchy- 
Riemann  conditions  are  valid.)  As  a  consequence,  the  final  solution  of  the  field  with  a 
radial  source  and  physical  blockage  is  perpendicular  to  the  field  with  an  angular  source 
and  physical  blockage,  resulting  in 

Next,  we  obtain 

^wr  _ _ 

~  Hr  ' 

Hence,  the  ratio  of  the  components  of  the  total  fields  from  the  angular  source  is  equal  to 
the  negative  of  the  ratio  of  the  fields  from  the  radial  source. 

The  fields  from  the  angular  source  were  easier  to  measure  in  the  laboratory.  The 
ratio  of  their  measured  components  was  checked  against  the  ratio  of  the  components  of 
the  fields  from  the  radial  sources,  which  is  easier  to  predict  theoretically,  as  the  model 
shows. 
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EXPERIMENTS 


The  experiment  was  set  up  to  measure  the  5  and  //  fields  from  a  wire  with  its  axis 
parallel  to  the  cylindrical  shell  blockage.  This  configuration  produced  an  angular  field, 
which  varied  as  -f\  ^om  the  source.  The  results  of  the  measurements  taken  at 
0  =  0  and  0  =  (figure  2)  show  a  negligible  amount  of  angular  field  disturbance,  which 
decreased  as  the  frequency  decreased.  The  model  predicts  zero  angular  disturbance 
effects  in  the  magnetostatic  limit.  Thus,  the  measurements  are  in  good  agreement  with 
the  model  for  this  experimental  configuration  and  these  angles.  The  slight  increase  in  the 
angular  field  disturbance  as  the  frequency  increased  can  be  explained  as  an  eddy  current- 
Lenz  law  response  of  the  physical  blockage  to  the  source  field. 

CONCLUSIONS  AND  RECOMMENDATIONS 

The  theoretical  model  for  the  angular  field  from  a  cable  source  predicted  no 
physical  distortion  of  the  angular  component  magnetic  field  at  0  =  0  and  d  =  n  with 
respect  to  the  axis  of  the  physical  blockage  (see  figure  1).  The  distortion  of  the  angular 
component  of  the  magnetic  field  was  calculated  by  taking  the  angular  derivative  of  the 
Legendre  polynomial  expansion.  After  these  derivatives  were  taken,  each  term  in  the 
Legendre  polynomial  expansion  had  a  sin  6  factor.  Hence,  the  disturbance  will  vanish  at 
0  =  0  and  6  =  %.  Experimental  measurements  of  the  angular  components  of  the 
magnetic  field  confirmed  this  behavior  (see  the  experimental  plots  in  figure  2). 

The  theoretical  model  also  predicted  a  periodicity  with  higher  harmonics  in  the 
disturbance  field,  and,  in  fact,  the  disturbance  field  will  vanish  altogether  on  certain  axes 
and  be  maximized  on  others.  Besides  the  angular  disturbance  caused  by  a  physical 
blockage  on  the  angular  field,  the  radial  disturbance  can  be  predicted  by  taking  the  radial 
derivatives  of  the  Legendre  polynomial  expansion. 

The  other  models  developed  were  for  the  vector  radial  field  somces  that  varied  in 
inverse  proportion  to  the  distance  and  cube  of  the  distance.  These  sources,  represented  in 
the  EMC  code  called  lEMCADS  as  power  supplies  and  transformers,  are  simplistically 
treated  by  the  lEMCADS  as  scalars;  however,  to  solve  for  the  distortion  effects  of  the 
physical  blockage  on  the  magnetic  source  fields,  the  sources  had  to  be  treated  as  vectors 
in  order  to  apply  the  boundary  conditions  that  require  the  angular  component  of  the  H 
field  and  the  normal  component  of  the  B  field  to  be  continuous  at  the  inner  and  outer 
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DISTANCE  (in.) 

Figure  2.  Measured  Magnetic  Field  From  a  Straight  Wire  Near  an  Iron  Pipe 


interfaces.  Upon  application  of  the  boundary  conditions,  a  set  of  10  equations  and  10 
unknowns  was  derived  and  solved  to  obtain  the  expansion  coefficients  for  the  Legendre 
polynomials,  which  resulted  in  the  total  field  demonstrated  in  this  report.  Although  these 
particular  models  were  not  checked  in  the  laboratory  since  the  mathematical  sources  are 
oversimplifications  of  real  power  supplies  and  transformers,  they  are  an  improvement 
over  the  mathematical  sources  used  in  lEMCADS  because  of  their  vector  (as  opposed  to 
scalar)  nature.  Realistically,  the  nearfields  from  transformers  and  power  supplies  have 
both  radial  and  angular  components  with  coefficients  in  spherical  coordinates.  The 
physical  blockage  for  the  models  assumes  that  the  angular  and  spherical  coordinates  of 
the  source  were  set  equal  to  it,  giving  a  worst-case  situation  when  compared  with  other 
coordinates. 


In  summary,  the  objective  of  this  investigation  was  to  develop  models  to  improve 
lEMCADS  by  incorporating  the  physical  blockage.  This  meant  keeping  the  source 
models  as  close  to  the  lEMCADS  models  as  possible.  We  could  not,  however,  use  the 
scalar  fields  designed  to  predict  worst-case  field  conditions  because  the  boundary 
conditions  at  the  physical  blockage  could  not  be  applied.  Although  the  complexity  of  the 
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nearfield  complicated  the  process,  we  reconstructed  sources  in  lEMCADS  to  give  them  a 
vector  nature  with  more  physical  reality. 


Future  work  should  focus  on  the  development  of  complete  nearfield  models  for 
transformers  and  power  supplies  rather  than  on  the  approximate  models  developed  here. 
The  same  techniques  described  earlier  for  determining  the  effects  of  physical  blockage 
should  be  applied  to  these  improved  source  models,  and  experiments  to  validate  the 
models  should  be  conducted.  The  source  model  used  for  the  cables  in  this  study  was 
exact  and  requires  no  improvement.  However,  experimental  measurements  of  the  effects 
of  the  physical  blockage  at  angles  others  than  0  =  0  and  6  =  %  should  be  undertaken. 
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Simplify  the  final  equations  for  the  coefficients  of  the  Legendre  Polynomials 
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Substitute  the  Coefficients 
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